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ABSTRACT 


The superconductor at the LaAIOa-SrTiOB interface provides a model system for the study of two-dimensional superconductivity 
in the dilute carrier density limit. Here we experimentally address the pairing mechanism in this superconductor. We extraot 
the electron-phonon spectral function from tunneling spectra and conclude, without ruling out contributions of further pairing 
channels, that electron-phonon mediated pairing is strong enough to aocount for the superoondueting critieal temperatures. 
Furthermore, we discuss the electron-phonon coupling in relation to the superconducting phase diagram. The electron-phonon 
spectral function is independent of the carrier density, except for a small part of the phase diagram in the underdoped region. 
The tunneling measurements reveal that the inorease of the chemical potential with increasing carrier density levels off and 
is zero in the overdoped region of the phase diagram. This indicates that the additionally induced carriers do not populate 
the band that hosts the superoondueting state and that the superoondueting order parameter therefore is weakened by the 
presence of charge carriers in another band. 


Introduction 

Interface superconductors are model systems for superconductivity in reduced dimensions' and provide input to the long¬ 
standing problem of understanding the mechanism of superconductivity in the layered cuprate superconductors. A prominent 
example is the superconductor at the LaAlOs-SrTiOs interface.In this two-dimensional electron liquid (2DEL), supercon¬ 
ductivity coexists with ferromagnetism'’"^ and spin-orbit coupling, providing the ingredients for exotic superconducting 
states such as hnite momentum pairing. '' Regarding the pairing mechanism, both conventional electron-phonon coupling 
and electronic pairing mechanisms are being considered.We recently mapped out the superconducting gap across the 
phase diagram and obtained a picture qualitatively similar to the phase diagram of the high-7^-cuprate superconductors: in 
the underdoped region the gap increases with charge carrier depletion.'^ This similarity between the doping dependence 
of the superconducting gap of LaAlOs-SrTiOa and the high-Tc-cuprates makes understanding the pairing mechanism in 
LaA 103 -SrTi 03 even more relevant. No experimental study of the pairing interaction exists, however. 

The LaA 103 -SrTi 03 interface 2DEL differs from the electron system in doped SrTi 03 because the band structures of the 
systems are different. But the critical temperatures are in the same range and it may well be that the pairing mechanism is the 
same. The superconductivity in doped SrTi 03 '® is generally explained in terms of the Bardeen Cooper Schrieffer (BCS) theory. 
Different phonon modes can contribute to the pairing interaction and the relative contributions of the different modes are still a 
topic of discussion.The critical temperature of doped SrTi 03 has a dome-shaped dependence on the carrier density n, 
with a maximum (400-600 mK) at an extremely small carrier density of ^lO'^ According to the BCS theory, the 

Tc of a superconductor depends on the Fermi surface area and on the electron-phonon coupling strength. As doped SrTi 03 
has a small Fermi surface area, it follows that the electron-phonon coupling would need to be exceptionally strong to explain 
superconductivity at such a small carrier density. Strong electron-phonon coupling is indeed possible because the Fermi energy 
and the plasma edge are smaller than some of the SrTi 03 phonon energies (predominantly the fourth longitudinal optical mode 
L04 at ^100 meV) involved. Therefore plasma excitations cannot screen these phonons well, yielding strong electron-phonon 
coupling, as observed in tunneling experiments.^^’At n > 5T0'^ cm"^, however, the plasma edge energy exceeds 100 meV, 
so that screening becomes more effective and the electron-phonon coupling is reduced. This mechanism has been proposed to 


explain the reduction of Tc of doped SrTi03 in the overdoped regime.^*' Furthermore, the SrTiOs longitudinal optical phonon 
modes have also been suggested to be important for the superconductivity in FeSe monolayers on SrTiOs.^^ 

To shed light into the pairing mechanism of the superconducting phase at the LaAlOs-SrTiOa interface we have performed 
tunnel experiments to spectroscopically measure the electron-phonon coupling a^F{(0) at the LaA103-SrTi03 interface. We 
observe coupling to the SrTi03 LO modes and conclude that these phonons likely provide the pairing interaction. We measured 
the evolution of both the chemical potential and the electron-phonon spectral function across the superconducting dome. The 
electron-phonon spectral function is found to only depend on the carrier density in the underdoped region. In the overdoped 
region the chemical potential is surprisingly constant. It is concluded that with increased doping the additionally induced charge 
carriers reside in a band that does not contribute to the pairing. Instead, the additional carriers result in Coulomb scattering of 
the electrons in the superconducting band, thereby reducing the superconducting gap. 

Results 

Figure 1 presents the differential conductance characteristic dl /dV (V) of a typical LaA103-SrTi03 tunnel junction at T = 
4.2 K. Here I is the tunnel current and V is the voltage applied between the 2DEL and the Au counterelectrode. The polarity 
of the voltage characterizes the sign of the interface voltage with respect to the top electrode bias; for V < 0 electrons tunnel 
out of the 2DEL. The tunnel characteristics are asymmetric, with a large tunnel conductance for V > 0 and a relatively small 
tunnel conductance for F < 0. At voltages well below the barrier height, the differential conductance of a tunnel junction 
is proportional to the density of states in the electrodes. Because the density of states of the Au electrode does not change 
significantly with energy, the d//dF (V) characteristic reflects the density of states of the 2DEL and the inelastic tunneling 
processes. Close to the conduction band edge of the 2DEL, the 2DEL density of states is strongly energy-dependent. When the 
absolute value of the negative voltage exceeds the chemical potential fi, the differential conductance is expected to almost 
vanish, as the tunnel conductance then probes the density of states in the bandgap. Here jj. is defined with respect to the band 
edge, so at y = 0 V the 2DEL density of states at £ = /r is probed. The differential conductance rapidly decreases for V < 0 
and a minimum is observed at V = —30 mV (Eig. 1). We identify this minimum with the conduction-band edge and attribute 
the increase of the conductance at larger negative bias voltages to inelastic tunneling processes and to barrier effects. In this 
part of the characteristic two prominent peaks are present, at V « —60 mV and V « —100 mV. 

To analyze the origin of the peaks in the inelastic tunneling conductance, the second derivative of the I{V) characteristics is 
analyzed, as the peaks in the second derivative correspond to the energies of the interacting boson modes.Eigure 2 presents 
the -d^//dy^(y) characteristics of seven tunnel junctions on different samples. Five of the samples were fabricated with standard 
LaA103-SrTi03 interfaces, while sample T6 was fabricated using a SrTi**03 substrate. The five standard samples have very 
similar characteristics, the peaks being at identical voltages. The four LO phonon modes of SrTi03 are at energies of 19.8, 
33.0, 58.8 and 98.6 meV (measured at room temperature).At the energies of the L02, L03 and L04 modes we indeed 
observe peaks in the -d^//dy^(y) characteristics. The LOl mode is not directly observed. Additional peaks are observed at, 
e.g., 77.9, 157.5, and 196.2 mV. These can be identified as harmonics of the phonon energies, L01-I-L03, L03-I-L04, and 
twice L04, respectively. The phonon energies extracted from the tunneling data are in good agreement with those observed by 
hyper-Raman measurements,^^ as listed in Table 1. The small difference originates presumably from the different measurement 
temperatures. One sample was grown using a SrTi^^03 single crystal.In this oxygen-isotope substituted sample, small shifts 
of the phonon energies can be expected. We observed a shift of 3.2 meV towards lower energy in the L04 mode. The other 
modes are shifted by less than 1 meV. The isotope effect is stronger for the L04 mode as this mode involves large displacements 
of the oxygen ions. Because similar shifts have been observed in Raman measurements on SrTi'^03 (not shown), a significant 
fraction of (larger than 60%) has to be present in the interfacial region of the sample, even after the LaA103 growth and 
annealing in a environment. 

Having shown that electron-phonon coupling can be directly observed in the LaA103-SrTi03 tunnel junctions, we now 
discuss its relation to superconductivity. We measured the gate-voltage dependence of the phonon-assisted tunneling in a device 
for which we previously determined the superconducting phase diagram.*^ A positive (negative) gate voltage accumulates 
carriers at (depletes carriers from) the 2DEL. Optimum doping with maximum is achieved at Vq = 0 V. At positive (negative) 
gate bias the system is overdoped (underdoped). The gate-voltage dependence of the tunnel characteristics is presented in 
Fig 3. For positive gate voltages the characteristics are not affected by Vq. Negative gate voltages, however, decrease the 
tunnel conductance significantly. This decrease is due to a change in the chemical potential ji that reduces the occupied density 
of states in the 2DEL. To precisely determine the gate-voltage dependence of jj., we analyzed the shifts in voltage in the 
conductance curves at positive voltages, see Fig. 4. These shifts are constant over a large voltage range and can be used to 
accurately determine A/r, the change in chemical potential in comparison to the overdoped cases. We determine /r for Vg > 50 
V by finding the crossing point between the almost constant conductance in the range -55 < V < -35 mV (mostly inelastic 
tunneling) and the strongly voltage dependence conductance in the range -25 < V < 0 mV (elastic tunneling). This yields /r = 
30 ± 2 mV, in good agreement with the data in Fig. 1 . The tunneling spectra in Fig. 3 show a small reduction of the conductance 
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close to the Fermi energy at V > -3 mV. In related cases, this rednction is attribnted to the Altshnler-Aronov correction to 
the density of states of an electron system with electron-electron interactions.^^ A detailed analysis of the Altshnler-Aronov 
correction is beyond the scope of this article. The most important observation is that the phonon-assisted tunneling peaks are 
consistently observed for all gate voltages. 

For quantitative analysis of the electron-phonon coupling, the inelastic tunneling probability has to be evaluated at the 
different energies. In the case of tunneling with the emission of real phonons studied here (as opposed to the case of the virtual 
phonon-coupling induced self-energy correction studied in superconducting tunnel junctions^^), the tunneling probability is 
proportional to The tunneling conductance is a function of both the density of states and the inelastic tunneling 

probability. The observed shape of the inelastic tunneling peaks is due to a convolution of the electron-phonon spectral function 
with the occupied density of states of the 2DEL. Because the voltage range in which elastic tunneling occurs (V > —pile) and 
the voltage range in which inelastic tunneling occurs (V < —file) are separated in energy, we have experimental access to the 
occupied density of states of the 2DEL. Here e is the electron charge. We therefore deconvoluted the density of occupied states 
(as measured by elastic tunneling in the voltage range -jije <V < 0) from the inelastic tunneling conductance (see Methods). 
This procedure yields a function proportional to the electron-phonon spectral function, which is shown in Fig. 5. The magnitude 
of the function has been normalized such that a^F (a))-da) reflects the ratio of the inelastic tunneling transmission in an energy 
range dm around co and the total elastic tunnel transmission. The dominant features of a^F (m) are the strong coupling at the 
L03 (~60 meV) and L04 (^100 meV) phonon modes. The phonon energies obtained from a^F (m), see Table 1, are are in 
good agreement with those obtained from the peaks in -d^//dV^(y). Next to the peaks from the phonons, a background that 
increases with increasing energy is present. In tnnneling, the barrier height decreases with increasing voltage and therefore 
some additional elastic tunneling is also present at voltages V < —file. The deconvolution procedure ignores this and the elastic 
part of the tunneling conductance results in the background. The a^F (m) function is unaffected by the gate voltage, except for 
an overall increase at negative gate voltages. 

Discussion 

The main objective of our study is to identify the pairing mechanism of the LaA103-SrTi03 2DEL by measuring the coupling 
of the electrons to bosonic modes with inelastic tunneling spectroscopy. Up to energies of 200 meV we hnd only coupling to 
the LO phonons of SrTi03. The phonon-assisted tunneling conductance in the 2DEL junctions is signihcantly larger than that 
observed in, e.g., Pb junctions^^"^® and is of similar magnitude as that observed in doped SrTi03 junctions.This indicates 
that the electron-phonon scattering crosssection is large. The measurements yield a^F (m) with an unknown proportionality 
constant and we therefore cannot calculate the critical temperature. So, the results are not nmambignons proof that the LO 
phonons do provide the pairing mechanism. Because we observe coupling to the LO phonons and do not observe coupling to 
other modes, we have to conclnde, however, that these tunnel spectroscopy measurements point clearly to electron-phonon 
coupling as the pairing channel. Note that next to the observed coupling to the LO phonons, coupling to the acoustic phonons 
may be present as well, as discussed for example in ref Our measurements namely do not discriminate coupling to acoustic 
phonons modes from the elastic tunneling conductance if the bosonic spectral function is not strongly energy-dependent. 

We next discuss the gate-voltage dependence of the electron-phonon coupling. The coupling strength is characterized by 
the McMillan parameter A that is obtained from a^F {(o) by^*^ 

X = 2f a^F{co)/co-dco. (1) 

Jo 

We extracted A from the data in Fig. 5 by integrating equation 1 over the energy range 30 < E < 145 meV. The resulting 
values for A are presented in Fig. 6a, where they have been normalized to those in the overdoped region. Approximately 60% 
of the conpling is dne to the L04 mode and approximately 25% of the conpling is dne to the L03 mode, independent of Vq. 
The doping dependence of A can be compared to that of the previously determined Tgap values, the temperatures at which the 
superconducting gap closes. A is constant for Vg > —50 V and increases with decreasing carrier density in the underdoped 
region of the phase diagram. In the underdoped region of the superconducting phase diagram, the increase of A with decreasing 
carrier concentration is qualitatively consistent with the increase of Tgap- However, in the main part of the phase diagram Tgap 
depends strongly on the applied gate voltage and A is constant. The decrease of Tgap in the optimally doped and overdoped 
region is therefore puzzling. In the following we show that this decrease can be directly nnderstood by considering the band 
structure of the 2DEL. 

Density functional theory calculations,®’^*’"*^ transport properties*^’and recent angle-resolved-photoemission (ARPES) 
data*^ indicate that the 2DEL comprises several bands: small electron mass d^^ orbital derived bands and large electron 
mass dxz and dy^ orbital derived bands. Becanse the dxz and dy^ orbital derived bands have larger momenta in the tunnel 
direction, tunneling occurs predominantly to those bands. This is consistent with the band structure determined from ARPES 
measurements:*^ the bottom of the d^j, and dyj orbital derived bands lies approximately 50 meV below the Fermi energy, in 
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reasonable agreement with the /r = 30 mV observed in tunneling. The bottom of the d-xy orbital derived bands lies approximately 
300 meV below the Fermi energy. Because the superconducting gap is observed in the tunneling characteristics with the 
temperature dependence of a primary order parameter,*^ we conclude that the and dy^ orbital derived bands host the 
dominant contribution to the superconductivity. 

Figure 6b presents the gate voltage dependence of the chemical potential in these bands. In the underdoped region jj. steadily 
increases with increasing carrier density, as expected. However, at optimal doping the increase of /r levels off and jj. is almost 
constant in the entire overdoped region. Because the shape of the tunneling characteristic is virtually independent of Vq, the 
additionally induced charge carriers have to reside in one or several higher-energy bands with large density of states that are not 
accessible to the tunneling, such as in a d^y orbital derived band or in a band further away from the interface. In agreement with 
this conclusion, signatures of such a band have been observed in transport studies,exactly appearing at gate voltages larger 
than the one at optimum (see also the Hall effect data in references'®’'*^). The charge carriers in this additional band are 
presumably not superconducting (or only superconducting due to the proximity effect from the other bands), but will provide 
more Coulomb scattering to the electrons in the superconducting bands, thereby explaining the reduction of superconductivity 
in the overdoped region. Electronic phase separation in superconducting and non-superconducting regions could also explain 
the doping independence of fi, but in this scenario the density of states is expected to change with doping, contrary to the 
measurements. 

In summary, we performed tunneling experiments to identify the superconducting pairing mechanism in the LaAlOs-SrTiOs 
2DEL. We determined a^F (co) and observed electron-phonon coupling across the entire superconducting phase diagram. 
We only observed coupling to the LO phonons of SrTi 03 and the coupling to the L04 mode was measured to be particularly 
pronounced. We conclude that electron-phonon coupling likely provides the pairing mechanism for the superconductivity. 
In the underdoped region the decrease of Tgap with increasing carrier density is possibly explained by a reduction of the 
electron-phonon coupling strength, but in the optimally doped and overdoped regions the electron-phonon coupling is doping- 
independent. In these regions an additional band becomes populated, as evidenced by tunneling measurements of the chemical 
potential in the 2DEL. The charge carriers in this band result in additional Coulomb scattering and thereby weaken the 
superconductivity, causing the reduction of Tgap. This scenario is intriguingly similar to the reduction of in overdoped 
gate-tuned MoS 2 ^^’^^ and to the constant chemical potential observed In cuprate Interface superconductors.®® 


Methods 

Experimental 

The tunnel junctions were fabricated by first growing a 4 or 5 unit cell thick layer of LaAlOs on Ti02 terminated®* SrTiOa by 
pulsed laser deposition to create the 2DEL.®^ Then a gold top electrode was deposited on the LaAlOs in situ. The gold layer 
was patterned using standard photolithography and selective chemical etching with a KJh-I 2 solution. In a final processing step 
ohmic contacts to the electron system were made by argon ion milling and Ti sputtering. The fabrication and characterization of 
the devices is described in more detail elsewhere.*® We fabricated devices using both SrTi *®03 and SrTi*^ 03 ®‘* single crystals. 
We did not see an effect of the oxygen isotope exchange on the normal-state transport properties or on the superconducting 
properties of the 2DEL. 

Extracting the electron-phonon spectral function 

In superconducting mnnel junctions, the self-energy correction to the BCS density of states due to the electron-phonon 
coupling can be directly observed and quantitatively modeled using the Eliashberg theory. In the LaA 103 -SrTi 03 2DEL this 
does not work because the self-energy correction is not observed. The phonon-assisted tunnel spectra are identical in the 
superconducting and normal state (not shown) because the superconducting gap (^50 /reV) is much smaller than the phonon 
energies involved. Therefore the normal state density of states has to be used to quantify the electron-phonon coupling. In 
normal metal phonon-assisted tunneling, the inelastic part of the tunneling conductance is a convolution of the electron-phonon 
spectral function and the density of states of the electron system. 

giiE) = kJ^ a^F{E - T)Aocc(-T)dT, (2) 

where gi is the inelastic conductance, is a constant and Nocc{E) is the occupied density of states. In case the density of states 
is constant in the energy range of interest, Nocc is the heaviside step function. Then the derivative of equation 2 with respect to 
the energy yields the proportionality between dgi{E)/dE and a'^F{(o). This relation is generally used in inelastic tunneling 
spectroscopy.^*’In case the density of states is not constant, changes in gi{E) are either due to changes in a^F{(o) or 
to changes in Nocc{E). When Nocc{E) is known, a^F{(o) can be extracted from gi{E) with the procedure described in the 
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following. A discrete version of Eq. 2 can be written as 


= AEitm ■ o) 

Here AE is the step size in energy. The summation is cut off after T = /r because the density of states is zero at larger energies. 
We identify the tunnel conductance in the range -/i/e <V < 0 as Nocc{E) and the tunnel conductance in the range V < -pLit as 
g-i{E). Now (/i + AE) can be used to determine (/i -t- AE ), because a?'F{E) is assumed to be zero for £ < /i. Following 
this, a^F{n -f 2A£) can be obtained. The procedure works best when the tunnel conductance at £ = /i is zero, giving a clear 
separation between density of states and inelastic tunneling. If this is not the case, a large spike will result in a^£(£) at £ = 
/i + A£. This spike can be removed by adjusting a constant value for a}F{E) for £ < jx. 
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Table 1 . The energies of the longitudinal optical phonon modes of SrTi 03 as observed by hyper-Raman scattering,and the 
present tunneling spectroscopy. The Raman data were taken at 300 K and the tunneling data were taken below 4.2 K. The 
numbers in parentheses indicate indirect measurements, derived from higher order processes. 


Phonon 

mode 

Hyper-Raman 

Tunneling 

-d2//dy2 

Tunneling 
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33.0 meV 
58.8 meV 
98.6 meV 
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60.4 meV 
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Figure 1 . d//dV (V) tunneling characteristic measured at T = 4.2 K. The dashed green line separates the voltage ranges where 
elastic and inelastic tunneling dominate in the transport. For F > 0 electrons tunnel from the Au electrode into unoccupied 
states in the 2DEL and for V <0 electrons tunnel from the 2DEL into the Au. In the inelastic tunneling regime, the transport is 
mainly by phonon-assisted tunneling with SrTiOs phonon modes. This process is illustrated in the density of states versus 
energy diagram of the junction in the inset, here the density of states of Au is shown on the left, that of the 2DEL on the right. 
If driven by voltages so large that their energy exceed phonon energies, electrons in the occupied states of the 2DEL (blue area) 
tunnel directly into the Au or they first make a transition to a virmal state in the gap by emitting a phonon (dashed blue area) 
and then transfer into the Au. Note that the linear form of the DOSIE) shown in the schematic is not the actual DOSIE) of the 
2DEL. 
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Figure 2. The second derivative -d^//dV^(y) characteristics of tunnel junctions showing phonon-assisted tunneling. The data 
were normalized to a maximum of 1 and offset by multiples of 0.5. Samples Tl, T2, T3, T4 and T6 have a 4 unit cell thick 
LaAlOs layer, while sample T5 has a five unit cell thick layer. The colored solid lines mark the energies of the LO phonon 
modes. Dashed lines mark sums and harmonics of these energies with the same color code. The two devices on 
LaA 103 -SrTi ^^03 sample are referred to by T6a and T6b. 



y(V) 


Figure 3. Gate-voltage dependence of the phonon-assisted tunneling. The tunneling characteristics dZ/dV (V) were measured 
at 0.05 K. 
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Figure 4. Gate voltage-dependent conductance characteristics at positive voltages, measured at 0.05 K. The shift in voltage in 
these characteristics was used to determine the change in chemical potential A/i that is shown in the inset. 



Figure 5. The electron-phonon spectral function a^F{(o) as calculated from the dZ/dV(V) data in Fig. 3. Each curve has 
been offset for clarity. 
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Figure 6. Gate voltage dependence of a) A, b) jj. and c) Tg^p and T^. The values for A are normalized to those in the overdoped 
region. In the underdoped region (brown background) the doping dependence of Tgap coincides with a change in A. In the 
optimally doped and overdoped regions (blue background) the reduction of Tgap coincides with a lack of increase of /r with 
increasing carrier density, indicating the additional charge carriers do not reside in the band that is superconducting, but 
interfere. The error margins of A reflect the uncertainty in the background subtraction procedure. The values for Tgap and Tc are 
taken from an earlier publication.'^ 
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